Integrative Biology

Glioblastoma (GBM) is a highly aggressive brain tumor with a five-year survival rate of less than 5%. The current standard of care established 20 years ago includes maximal surgical resection and administration of alkylating agent temozolomide (TMZ). GBM is highly invasive, and GBM cells that evade surgical resection can become resistant to TMZ and develop new aggressive secondary tumors. Post-relapse there are few treatment options available to patients. Tissue engineering approaches suggest the opportunity to develop in vitro models of the GBM tumor microenvironment that may accelerate the discovery of novel therapies for GBM. Here, we report the adaptation of hydrogel microdroplets (microgels) to encapsulate GBM cells in a tailorable 3D matrix to assess patterns of growth and to screen TMZ drug response using patient-derived xenograft (PDX) specimens. We exploit a unique aspect of the microgel system to account for the cellular heterogeneity within the tumor microenvironment (TME). We combine cell-laden microgels generated from TMZ-resistant and TMZ responsive variants of the same PDX specimens to create heterogeneous populations with varying levels of drug sensitivity. We demonstrate a range of drug resistance phenotypes as a function of the ratio of TMZ-responsive to resistance cells and identify the population required for TMZ-resistance to overtake take the response. We then investigate the influence of tumor mimetic shifts in hyaluronic acid bioavailability and hypoxia on patterns of TMZ resistance. We show exposure to matrix-bound hyaluronan increases TMZ resistance and the glioma stem cell population in both cell variants. Lastly, we report an increase in TMZ sensitivity but divergent changes in the GSC subfraction for TMZ resistant vs responsive GBM in the presence of hypoxia. Together, we demonstrate the versatility of cell-laden microgel approach to replicate heterogenous tumor populations, model shifts in the tumor microenvironment, and rapidly screen therapeutic response.

High-risk neuroblastoma patients face poor survival despite intensive treatment. Drug testing using patient-derived models can support therapy prioritization for precision medicine and drug development. Models incorporating tumour microenvironmental components, such as co-cultures and short-term cultured patient samples containing substantial non-malignant cell fractions, could better recapitulate microenvironment-dependent drug responses. However, conventional viability assays measure the combined signal from all viable cells in a well and therefore cannot determine tumour-specific drug responses. Here, we establish a microscopy-based readout to quantify cell-type-specific viability in two complementary settings: neuroblastoma-PBMC co-cultures and freshly dissociated patient tumour samples. In the co-cultures, PBMCs were pre-labelled with a cell-tracking dye, and Calcein staining was used to independently quantify the viability of tumour cells and PBMCs in the same well. The Calcein-based viability readout correlated strongly with conventional CellTiter-Glo measurements and was compatible with automated high-throughput drug screening. The imaging workflow enabled identification of compounds with differential efficacy in co-culture versus monoculture and distinguished tumour-specific effects from PBMC toxicity. The microscopy-based viability readout was further adapted to short-term cultured patient samples. Neuroblastoma tumour cells were distinguished from the non-malignant cells using a combination of tumour-specific surface markers NCAM, L1CAM and B7H3. This enabled determination of tumour fractions and measurement of tumour-specific drug responses. Tumour fractions varied substantially between patient samples, highlighting the importance of tumour-specific viability measurements. Together, the microscopy-based viability readout for co-cultures and patient samples enables scalable assessment of tumour-specific drug responses.

The endothelial cells (ECs) that line blood vessels continuously sense and respond to the physical forces exerted by blood flow. In vivo, pulsatile arterial flow interacts with vessel curvature, branching and other anatomical features to generate complex local hemodynamic environments that dictate the magnitude, direction, pulsatility, and oscillatory nature of wall shear stress experienced by ECs. Currently, accessible and reproducible in vitro models of complex pulsatile flow that recapitulate in vivo vascular anatomy remain limited. Here, we combine a novel rotational-flow endothelial culture platform with detailed computational fluid dynamics (CFD) modeling to characterize four well geometries designed to generate distinct hemodynamic environments. CFD analyses demonstrate that these geometries intrinsically generate pulsatile flow and produce reproducible spatially distinct regions of wall shear stress magnitude, pulsatility, and oscillatory shear within a single culture well. Endothelial alignment mapping and functional assays reveal region-specific cellular responses to the predicted local flow conditions that closely corresponded to the predicted local hemodynamic environment, linking complex flow patterns to endothelial adaptation. The technical advancements of our modeling efforts should support a faster, cheaper, simpler, and--importantly--validated framework for future investigation into EC mechanobiology under complex flow conditions. HIGHLIGHTSO_LISimple engineered well geometries generate distinct hemodynamic microenvironments, mimicking in vivo vascular structures, using a conventional orbital shaker. C_LIO_LIComputational fluid dynamics (CFD) reveals spatially distinct patterns of wall shear stress, pulsatility, and oscillatory shear applied to ECs within individual culture wells. C_LIO_LIHigh average wall shear stress and elevated oscillatory shear index induces a unique perpendicular alignment of ECs to the dominant flow vector. C_LI

Control of fluids is a hallmark of microfluidic systems and fundamental for the successful application of microfluidic devices. Trigger valves use geometric features to autonomously control the release of fluids in microfluidic devices. Our previous work has adapted geometries used in closed trigger valve systems to enable use in open systems, allowing for open microfluidic devices with up to three trigger valves. Here, we focus on the parallel co-flows produced by sequential release of trigger valves and present a model that predicts their layer widths as a function of the geometric characteristics of the different side channels of each trigger valve. We show layered co-flows with widths as low as 50 microns. Additionally, we expand the use of trigger valves in open microfluidic devices by incorporating 1) varied step heights, 2) devices with up to seven trigger valves, and 3) use of varied fluids and plastics. To validate the implementation and use of these trigger valves in open systems, we have developed a theoretical framework to compare predicted outcomes (i.e., fluid travel distance, velocity, and layering width) with our experimental values. This theoretical work offers applications in various fields, including hydrogel patterning for 3D cell culture, organ-on-a-chip models, at-home sample preparation, and autonomous microfluidic systems for biosensing.

Intracellular iron homeostasis is increasingly recognized as a regulator of cancer cell behavior, but how iron distribution influences extracellular matrix (ECM) organization and invasion remains poorly understood. Here, we show that loss of divalent metal transporter 1 (DMT1/SLC11A2) disrupts intracellular iron homeostasis and promotes cancer cell invasion through an iron-ER-ECM axis. In MDA-MB-231 cells, DMT1 knockout (KO) reduced total iron content but increased the labile iron pool (LIP) in both 2D and 3D culture models, indicating altered intracellular iron distribution. Across transcriptomic and phenotypic readouts, DMT1-dependent effects were more evident in 3D than in 2D models, with DMT1 KO inducing endoplasmic reticulum (ER) stress and impaired collagen/ECM organization. Functionally, the DMT1-loss phenotype was marked by reduced 2D motility, whereas in 3D spheroid models DMT1 KO cells displayed enhanced invasive outgrowth in both Matrigel and collagen matrices. Iron chelation further modulated this phenotype in a DMT1-dependent manner. Pharmacologic induction of ER stress phenocopied the loose spheroid architecture and invasive behavior, supporting ER stress as a mechanistic link between altered iron handling and ECM destabilization. Together, these findings identify intracellular iron distribution, rather than total iron abundance alone, as a determinant of ECM integrity and context-dependent cancer cell invasion. Significance StatementOur study identifies an iron-ER-ECM axis through which intracellular iron homeostasis regulates cancer cell invasion. Total cellular iron content alone is insufficient to predict invasive behavior without considering how iron is distributed within the cell. By preserving intracellular iron homeostasis and ER function, DMT1 supports collagen synthesis and maintains ECM integrity. In contrast, DMT1 loss disrupts these processes, promoting formation of loosely aggregated spheroids and enhanced invasion in 3D tumor models despite reduced total iron levels. These findings challenge the assumption that lowering bulk iron uniformly suppresses invasive phenotypes and instead highlight intracellular iron trafficking as a potential therapeutic target for limiting cancer cell invasion.

Integrin v{beta}8 is a major activator of latent transforming growth factor-{beta} (TGF-{beta}) and an emerging therapeutic target in cancer and immune regulation. However, functional assays that directly measure v{beta}8-mediated activation of latent TGF-{beta} in a physiologically relevant context remain limited. Here, we report a co-culture cell-based reporter assay for quantitative measurement of v{beta}8-mediated activation of latent TGF-{beta}1. NIH/3T3 reporter cells were engineered to express a SMAD-responsive NanoLuc reporter, constitutive firefly luciferase for internal normalization, and cell-surface GARP-latent TGF-{beta}1. When co-cultured with v{beta}8-expressing LN-229 cells, reporter cells produced a robust signal that directly reflected localized latent TGF-{beta}1 activation. The assay demonstrated stable expression of the required biological components, reproducible signal-to-background performance, and sensitivity to benchmark v{beta}8-blocking antibodies. Inhibition studies showed potent dose-dependent blockade by an anti-v{beta}8 antibody. In contrast, pan-TGF-{beta} neutralizing antibody displayed markedly weaker apparent potency, suggesting that targeting localized v{beta}8-mediated activation is more effective than neutralizing released TGF-{beta} in this assay context. The assay also enabled screening and ranking of anti-v{beta}8 antibodies, identifying several high-potency clones, and detected v{beta}8-mediated activation of a non-cleavable latent TGF-{beta}1 mutant. This platform provides a sensitive, internally normalized, and scalable approach for mechanistic studies and therapeutic discovery targeting the v{beta}8-TGF-{beta} axis.

AO_SCPLOWBSTRACTC_SCPLOWCell fate decisions are driven by gene regulatory networks (GRNs). While the mutually inhibitory toggle switch effectively models binary fate decisions, fully connected inhibitory networks with more than two nodes fail to capture multi-fate decisions due to the low prevalence of "single high states", where only a single master regulator is highly expressed. The goal of this study is to find network structures that support all single high states. We find that the only network that attains the highest possible prevalence of all single high states within the set of monotone Boolean (MB) models is completely disconnected. Since biological networks typically require connectivity, we investigate network structures that support equipotency, where all single high states have equal prevalence within MB models. Finally, we characterize the networks that support multistability between all single high states, finding that it is possible only in networks in which each node either has self-activations or is inhibited by every other network node. Our findings provide a theoretical framework for understanding the network design principles that can support simultaneous differentiation into multiple distinct cell types.

High-grade serous carcinoma (HGSOC) is the deadliest subtype of ovarian cancer, characterized by high metastatic rates. HGSOC is typically diagnosed at late stages, and treatment options are limited, resulting in a 60% recurrence rate. HGSOC cells exhibit metabolic plasticity, dynamically shifting between glycolysis and oxidative phosphorylation (OXPHOS) to meet energy demands for tumor progression. To evaluate therapeutic strategies that target metabolic vulnerabilities, we developed a microphysiological system (MPS) that recapitulates the heterogenous cell states and bioenergetic distribution of HGSOC solid tumors. Our platform utilized HGSOC spheroids embedded in a collagen hydrogel that mimics the extracellular matrix to capture tumor progression in the ovary. We used atovaquone (ATO), an FDA-approved OXPHOS inhibitor, to prototype the capabilities of our platform to investigate metabolic plasticity in HGSOC. Treatment with ATO decreased viability and invasion of HGSOC spheroids. Crucially, ATO exhibited no cytotoxicity toward biomimetic blood vessels, preserving their integrity and permeability. Metabolic imaging revealed that ATO induces an oxidative state in the outer region of the spheroids. At the invasive front, ATO disrupted mitochondrial organization, forcing collective cell migration and eventually inducing breakdown of mitochondrial networks. Furthermore, ATO decreased YAP/TAZ pathway activity in the outer region of the spheroid, providing a potential mechanism for hindered cell invasion. Collectively, our data demonstrates that a low-potency OXPHOS inhibitor like ATO can effectively target metabolic plasticity to suppress HGSOC spheroid progression. Overall, this platform successfully recapitulated metabolic heterogeneity and provided a workflow for safely testing other drugs that target cancer metabolism.

The suspension culture of pluripotent stem (PS) cells in stirred bioreactors poses a delicate balance between maintaining homogeneous cell dispersion and avoiding excessive shear stress that can compromise cell viability and pluripotency. In this study, we used computational fluid dynamics (CFD) coupled with a discrete particle method (DPM) to simulate iPS cell behavior in a 5 mL delta-impeller stirred tank. Our analysis revealed that upward flow at the tank bottom and downward flow at the top are critical for maintaining a stable suspension. To optimize the stirring protocol, we applied Bayesian optimization to identify a time-dependent stirring schedule that begins with a high-speed phase for resuspension, followed by a low-speed phase for sustained suspension with minimal hydrodynamic stress. The optimized schedule demonstrated improved suspension ratio and reduced slip velocity, indicating lower mechanical stress on cells. These findings provide engineering insights into scalable bioreactor operation, contributing to the design of robust iPS cell manufacturing systems.

Adoptive T-cell therapies rely on the identification and expansion of rare tumor-reactive T cells, yet current enrichment strategies are limited by the low abundance of these cells and complexity of their functional enrichment. Here, we present a microfluidic platform that exploits hydrodynamic shear as a controllable parameter for enriching antigen-specific T cells through peptide-major histocompatibility complex (pMHC)-mediated capture. An eight-channel microfluidic device was engineered to simultaneously interrogate a range of wall shear stresses while maintaining uniform cell delivery, enabling systematic identification of shear conditions that maximize antigen-specific enrichment. Using engineered MART-1-specific Jurkat cells, we demonstrate that T-cell capture is jointly regulated by wall shear stress and pMHC density, with intermediate shear preferentially enriching antigen-specific cells over nonspecific binders. Translation of the optimal operating condition to a high-throughput single-shear device enabled approximately 35-fold enrichment of antigen-specific T cells from peripheral blood mononuclear cells containing only 0.05% target cells. We further show that peptide-MHC complexes isolated directly from melanoma whole-cell lysates support shear-dependent enrichment comparable to recombinant pMHCs. Finally, primary MART-1-specific CD8 T cells enriched using tumor-derived pMHCs retained the ability to recognize melanoma cells and upregulated the activation marker CD137 following antigen-specific stimulation. Together, these findings establish hydrodynamic shear as an orthogonal parameter for antigen-specific T-cell enrichment and provide a framework for integrating force-based selection with tumor-derived pMHCs to isolate functional antigen-specific T cells using tumor-derived pMHCs.

Mesenchymal Stromal Cells (MSC) face several heterogeneity challenges hindering clinical and commercial success. Employing a multiple response model, interplay between donor heterogeneity, and critical processing parameters (CPPs), effects on MSC potency and cell expansion attributes were investigated through computed composite attribute scores. Twelve unique CPP combinations were tested in thirteen marrow-derived MSC(M) and five adipose-tissue MSC(AT) training and test datasets, respectively. Donor heterogeneity and select CPP conditions affected a curated gene panel (surrogate for MSC potency); while MSC expansion was primarily influenced by CPPs. Model performances were evaluated against clinical effectiveness data from a previously deployed clinical trial; top-performing model predicted donor rankings coincided with clinical effectiveness data, validating the modeling approach used. Our model predicted that only 8% of tested donors were agnostic to CPPs; a majority (62%) of donors showed CPP-dependent optimal composite quality attributes, with MSC seeding density as a key driver; medium supplementation and oxygen preferences were highly donor dependent. Approximately 30% of donors performed poorly at all conditions tested and may be prospectively identified using a subset of genes (TGFB, VEGF, PDCD1LG1, PDCD1LG2, IDO). Model predicted optimal parameters worked for 69% of tested donors, while sub-optimal parameters worked for only 23% of donors and were confirmed in an independent CD14+ macrophage assay. Our integrated computational and experimental framework predictably identified interactive effects of donor heterogeneity and CPP conditions to optimize MSC potency attributes.

Osteoarthritis (OA) drug development remains constrained by preclinical models that fail to recapitulate the multicellular interactions that regulate human joint inflammation and extracellular matrix degeneration in response to investigational drugs. Tanezumab, a humanized anti-nerve growth factor monoclonal antibody developed for non-opioid pain relief, advanced to late-stage clinical trials but was discontinued due to unresolved joint-localized safety concerns, including rapidly progressive OA. This study evaluated whether a human microfluidic joint-on-a-chip co-culture system could detect early biomarker responses to tanezumab exposure that were not apparent in conventional chondrocyte monoculture. Tanezumab was first tested in human chondrocyte monoculture under untreated and disease-like (i.e., IL-1{beta}-treated) conditions. Across a 20-analyte panel of inflammatory and matrix-remodeling biomarkers, statistically significant monoculture responses to tanezumab were limited to decreased IL-1{beta} from 335 to 132 pg/mL ([~]0.39-fold) and increased IL-8 from 575 to 675 pg/mL ([~]1.17-fold). Major OA-associated matrix-remodeling markers, including MMP-1, MMP-3, and MMP-13, remained largely unchanged, indicating that monoculture conditions are insufficiently sensitive to detect clinically predictive drug-related molecular changes. Tanezumab was then evaluated in co-cultures containing chondrocytes, osteoblasts, fibroblast-like cells, and macrophages under low-inflammation (i.e., M0 macrophage-based) and high-inflammation (i.e., M1 macrophage-based) conditions. In the M0-based co-culture, tanezumab increased MMP-1 from [~]4.20 x 104 to [~]6.20 x 104 pg/mL ([~]1.48-fold), MMP-3 from [~]8.00 x 104 to [~]1.20 x 105 pg/mL ([~]1.50-fold), and MCP-1 from 2.85 x 103 to 4.31 x 103 pg/mL ([~]1.51-fold). In contrast, the M1-based co-culture showed decreases in MMP-13 from [~]1.66 x 104 to [~]1.17 x 104 pg/mL ([~]0.70-fold) and IFN-{gamma} from [~]1.95 x 104 to [~]1.56 x 104 pg/mL ([~]0.80-fold), changes that may appear beneficial despite the drugs known clinical risks. Collectively, these findings show that low-inflammation multicellular co-culture revealed coordinated matrix remodeling and inflammatory responses to NGF blockade that were missed in monoculture and were partly obscured in highly stimulated disease-like conditions. This platform may provide a useful, human-relevant approach for safety signal assessment and early evaluation of OA therapeutics within a defined context of use focused on joint-specific, tissue-level drug-response testing.

Purpose: Hyperthermia is emerging as an adjunct strategy in pediatric oncology, yet its translation is limited by poor understanding of how different modalities impact complex tumor microenvironments. Neuroblastoma (NB), the most common extracranial solid tumor in children, displays profound heterogeneity that hampers therapeutic predictability. Here, we performed the first systematic head-to-head comparison of photothermal therapy (PTT) and magnetic hyperthermia (MH) in tissue-engineered NB (TE-NB) models. Methods: TE-NB scaffolds incorporating five NB cell lines were loaded with magnetic nanoparticles (MNPs) and subjected to PTT (808-nm laser, 130 W/cm2, 10 min) or MH (285 kHz, 20 mT, 60 min). Constructs were analyzed at 24 h, 48 h, and 5 d post-treatment for DNA content, cell viability, proliferation (Ki67 immunohistochemistry), and apoptosis (caspase-3/7 staining). Results: MH produced consistent MNP-dependent heating with minimal background, while PTT was dominated by nonspecific medium absorption. Both modalities modulated proliferation within 24 h, but effects varied sharply by cell line and scaffold region, reflecting microenvironmental heterogeneity. By 48 h, PTT often triggered paradoxical increases in proliferation, whereas MH disrupted scaffold integrity, reduced DNA content, and suppressed Ki67 expression. Notably, neither modality induced sustained caspase-3/7 activation, indicating that cytotoxicity proceeds via non-apoptotic pathways. Conclusion: Our findings position MH as a superior modality for uniform heating and proliferation control in 3D NB models, but also highlight that hyperthermia should be considered a context-dependent modulator rather than a binary cytotoxic agent. By integrating patient-specific TE-NB platforms, precision hyperthermia could define individualized therapeutic windows, optimize combinations with pro-apoptotic or immunomodulatory agents, and accelerate translation of hyperthermia strategies for children with NB.

Generating bioengineered cartilage that recapitulates the depth-dependent phenotype, structure, and function of native articular cartilage remains a challenge. While cartilage is rich in aggrecan and type II collagen, proper function depends on depth-dependent protein expression. Superficial zone chondrocytes (SZCs) secrete proteoglycan-4 (PRG4) to lubricate the cartilage surface. Deep zone chondrocytes produce type X collagen (COLX) to support compressive loading and load transfer to subchondral bone. We previously demonstrated that passaged full-thickness chondrocytes (FTCs) and zonal chondrocytes can re-express cartilage and zone-specific markers following scaffold-free three-dimensional (3D) culture in redifferentiation media. However, in the absence of an instructive matrix, cells expressed low levels of zone-specific proteins and exhibited limited depth-dependent organization. We hypothesize that synthetic extracellular matrix with zone-specific microenvironmental cues will guide zonal differentiation. To this end, passaged primary bovine chondrocytes were encapsulated in a soft, hyaluronan (HA)-based, cell-adhesive, and protease-degradable hydrogel established via bioorthogonal tetrazine (Tz) ligation with norbornene (Nb). When supplemented with TGF{beta}3, FTCs deposited aggrecan and type II collagen with minimal type I collagen. Application of interfacial tetrazine ligation with trans-cyclooctene (TCO) during cell culture resulted in matrix stiffening, leading to upregulation of COLX expression. Conversely, SZCs cultured in soft hydrogels exhibited the greatest PRG4 expression. Establishment of a trilayered construct with region-specific stiffness via the diffusion-controlled reaction promoted PRG4 and COLX expression in defined zones. Together, these findings demonstrate that tunable HA-based hydrogels can enhance zone-specific chondrocyte phenotypes and promote the formation of zonally organized cartilage.

Cell surface receptors are activated by ligand binding and transmit signals into the cell. Epidermal growth factor (EGF) receptor (EGFR) signaling regulates cell growth, differentiation, and survival, and its dysregulation is linked to cancer. Recent experiments show that ligand binding to EGFR is enhanced for receptors in tetraspanin nanodomains on the cell surface. We use kinetic modeling of receptor confinement, ligand binding, and internalization to compare confinement and signaling behavior for EGFR with spatially localized ligand binding to a hypothetical receptor that has uniform ligand binding anywhere on the cell surface. We find that introducing a membrane domain that confines and enhances ligand binding to receptors leads to more consistent confinement across ligand levels, raises necessary ligand levels for steady-state signaling, and flattens and extends the signaling response to sudden ligand concentration increases. This confining domain that enhances ligand binding provides the cell with a distinct regulatory mechanism to tune its signaling response. We also find that the concentration of receptors in signaling states and the fraction of receptors in signaling states respond to ligand at different ligand concentrations, with substantial increase of the concentration of receptors in signaling states occurring at a much lower ligand concentration than a substantial increase of the fraction of surface receptors in signaling states. This quantitative modeling of spatially restricted receptor activation applies to other receptors with similar characteristics and builds towards physical principles of receptor signaling.

Vascular endothelial cells (ECs) form a monolayer lining blood vessels and serve as a barrier between blood and tissues. Clathrin-mediated endocytosis (CME) is a major internalization pathway that involves a physical conformational change of the plasma membrane to form a vesicle and is therefore sensitive to the local environment. ECs are subjected to a myriad of fluid shear stress (FSS) rates from circulating blood, which we hypothesize affects CME. To test this, we used simultaneous two-wavelength axial ratiometry (STAR) microscopy, which provides nanoscale axial resolution, to determine the frequency and morphology of clathrin-coated vesicles as they form. Human umbilical vein endothelial cells (HUVECs) were transfected with dual-tagged clathrin light chain a (CLCa-iRFP-EGFP) and cultured under 10 dyn/cm2 FSS. CME activity was elevated in cells cultured under flow and assayed in static or flow conditions compared to statically cultured and imaged controls, indicating that FSS-induced changes to CME were maintained shortly after flow cessation. Single vesicle analysis showed cells cultured in FSS had a slight preference for vesicle formation with a flat-to-curved clathrin transition compared to control. Next, to assess the impact of different FSS rates, we cultured HUVECs at 20 and 40 dyn/cm2 FSS. We found total CME frequency was elevated compared to control at 20 dyn/cm2, but not 40 dyn/cm2. HUVECs cultured at both 20 and 40 dyn/cm2 had vesicles with increased lifetime and enhanced stability, as well as a higher proportion of vesicles formed through a flat-to-curved transition of clathrin.

Producing bone and cartilage in a controlled and localized manner remains a significant challenge in regenerative medicine. This study investigated the ability of keratin- and polyethylene glycol (PEG)-based degradable hydrogels to deliver bone morphogenetic protein 2 (BMP-2) and leukocyte cell-derived chemotaxin 1 (LECT-1; also known as chondromodulin-1) intramuscularly to induce ectopic tissue formation. Adult male CD-1 mice received intramuscular implants of keratin-PEG gels containing a fixed dose of BMP-2 and increasing amounts of LECT-1. After two weeks, implants and surrounding muscle were analyzed using computed tomography (CT) and histology. The results showed that BMP-2 is necessary for forming new bone and cartilage, whereas LECT-1 alone appeared to trigger muscle dedifferentiation without ossification or chondrogenesis. Co-delivery of BMP-2 and LECT-1 enhanced bone and cartilage formation in a dose-dependent manner: higher LECT-1 doses led to proportionally more ectopic cartilage (linear correlation, r2 {approx} 90%), while bone formation peaked at the third LECT-1 dose at approximately twice the volume of the BMP-2-only group. These findings indicate that muscle-resident cells may be capable of reverting and switching to mesenchymal lineages, recapitulating endochondral ossification. The platform offers a promising strategy for growing bone and cartilage autografts within skeletal muscle bundles.

Mechanistic ordinary differential equation models are widely used in systems biology to represent biochemical networks, population dynamics, cell-state transitions, and other biological processes; however, their predictive value depends critically on accurate parameter estimation from noisy and often sparse experimental data. In this tutorial, we present the Weak-form Estimation of Nonlinear Dynamics (WENDy) method as a forward-solver-free approach that reformulates parameter estimation as a covariance-corrected weak-form regression problem by integrating the model equations against compactly supported test functions. We present the background on the methodology through the lens of the familiar logistic equation, and we demonstrate applications of the method on real experimental data through two systems biology examples: a glycolytic oscillator with relatively dense time-course data and a sparse epithelial-mesenchymal cellstate transition model with multiple experimental replicates. Ultimately, using WENDy, we estimate interpretable biological parameters with uncertainty for systems with noisy and sometimes sparse available experimental data.

Cell states are increasingly conceptualized as attractors of high-dimensional dynamical systems, yet quantitative approaches for integrating phenotypic information into this framework remain limited. Here, we take an image-based approach that combines unsupervised machine learning (ML) with timelapse imaging to extract and characterize the temporal dynamics of morphological features. Using a cell line with endogenously tagged VE-cadherin, we acquired brightfield and fluorescence timelapse images of human induced pluripotent stem cell-derived endothelial cell (hiPSC-EC) monolayers, which adopt distinct phenotypes at two different magnitudes of shear stress in terms of their morphology, behavior, and VE-cadherin organization. To quantify these phenotypic cell states without segmentation, we trained a diffusion autoencoder to predict VE-cadherin signal from brightfield images. We identified interpretable ML-based features representing cell orientation, elongation, and local density. Treating these variables as dimensions of a morphological state space, we estimated a data-driven vector field and found that the two observed phenotypic cell states correspond to stable fixed points of the inferred dynamical system. Mapping measured cell migration coherence onto this space further distinguished the states. Imaging cells across intermediate shear stresses revealed a regime of bistability in which both states coexist, indicating that the shear-stress-dependent transition between endothelial cell states occurs as a bifurcation of the inferred dynamical system. Finally, we applied this method to study an N-terminal truncation of VE-cadherin, finding that mutated cells preserve alignment and coherent migration, but exhibit altered morphology and increased migration speed. This work demonstrates the applicability of a dynamical systems approach to quantitatively characterize morphological aspects of cell state from interpretable ML-based features.

Control over network dynamics at different length scales is a feature of natural materials challenging to replicate in synthetic hydrogels. Hydrogel viscoelasticity is commonly controlled by tuning the kinetics of reversible crosslinks; however, this strategy inherently links the resulting macroscale and nanoscale dynamics of the individual network components. Taking inspiration from biological materials that feature lipids as structural elements, we introduce Lipid Network Crosslinked (LINC) hydrogels that exploit the mobility of individual lipids within self-assembled liposomes as covalent, network-crosslinking points. These mobile, covalent crosslinks increase hydrogel stress relaxation rates over 20-fold compared to polymer-only hydrogels with equivalent crosslinking chemistries and stiffnesses. We demonstrate that liposome design parameters, including degree of surface functionalization and tail saturation, provide a means to independently control the macroscale storage moduli and stress relaxation behavior. Finally, as an application where control over network dynamics at different length scales is critical, we placed cell-adhesive ligands onto more mobile or less mobile network elements. Human neural progenitor cells cultured within LINC hydrogels of identical macroscale viscoelasticity significantly altered their phenotype in response to nanoscale ligand dynamics. These results establish LINC hydrogels as biomimetic materials that leverage nanoscale lipid mobility within a macroscale polymeric network to control dynamics at multiple length scales.